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ABSTRACT 

The thermal stability of the (020)2 cluster dimer consisting of two C20 fullerenes is exam- 
ined using a tight-binding approach. Molecular dynamics simulations of the (020)2 dimer at 
temperatures T = 2000 - 3500 K show that the finite lifetime r of this metastable system is 
determined by two fundamentally different processes, the decay of one of the C20 fullerenes 
and the fusion of two C20 fullerenes into the C40 cluster. The activation energies for these 
processes E a w 3.4 and 2.7 eV, respectively, as well as their frequency factors, have been 
determined by analyzing the dependence of r on T. 
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In 2000, after the experimental detection of the smallest fullerene C 2 o [1], Fig. la, on the 
surface of which C-C bonds form only pentagons, whereas hexagons are absent (in contrast to 
Cjv fullerenes with N > 20), a problem of synthesizing the condensed phase (fullerite) based 
on this fullerene by analogue with the known fullerite based on Ceo fullerenes [2, 3] arised. 
Interest in this problem is stimulated, because the C20 fullerite may be a high-temperature 
superconductor [4-6]. Indeed, preliminary calculations point to a possibility of formation of 
three-dimensional crystals consisting of C20 fullerenes [4-6]. There is yet disagreement on 
the crystal structure of such a hypothetical cluster matter [5,6], but all authors agree that 
C 2 o clusters in it should be connected with each other by strong covalent bonds. This is a 
fundamental difference of the C20 fullerite from the Cqq fullerite in which clusters are bonded 
with each other due to a weak van der Waals attraction. 

The strong intercluster interaction can significantly distort the shape of the C 2 o clusters in 
the fullerite with respect to the spherical shape of the C 2 o fullerene and, as a consequence, 
result in the transition to another crystalline (or even amorphous) modification. It seems to 
be reasonable to begin the analysis of the stability of the C 2 o fullerite with the analysis of 
the stability of the simplest system consisting of two interacting C20 fullerenes, the (020)2 
cluster dimer. Experimental indications of the existence of charged (020)2" dimers [as well as 
complexes (C2o)fcWith k — 3 - 13] were obtained in [7]. The dimerization of neutral fullerenes 
C 2 o was theoretically studied in [4] and, most comprehensively, in [8], where eight different 
(C 2 o)2 isomers were found, their structural and energy characteristics were calculated, and it 
was shown that the so-called open- [2+2] isomer has the lowest energy, see Fig. lb (numbers 
in the square brackets are the numbers of the atoms of each fullerene that are involved in the 
formation of intercluster bonds). 
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The aim of this work is to theoretically analyze the stability of the (C 2 o) 2 dimer with respect 
to its transition to another atomic configuration. It is known (see, e.g., [9]) that the specific 
(per atom) binding energy of carbon clusters, which is defined as the difference between the 
energy of an isolated carbon atom and the specific potential energy of the iV-atomic cluster, 
Eb(N)/N = E(l) — E pot (N)/N, increases with N (this is the property of the overwhelming 
majority of the elements of the Periodic Table, except for metastable clusters of nitrogen, 
helium, and some other elements [10, 11]). For this reason, it is energetically favorable 
for two (or several) carbon clusters to merge into one larger cluster, thereby increasing the 
binding energy [10] (i.e., reducing the potential energy), by analogy with fusion of light 
atomic nuclei. At the same time, our previous analysis [12] of the thermal stability of the 
isolated C 2 o fullerene showed that, when heated, it passes to various energetically unfavorable 
configurations with lower E^ values rather than to a more stable isomer with the maximum 
binding energy ("bowl"). This is due to the ability of carbon structures to form numerous 
intermediate metastable states separated by high energy barriers from configurations with a 
lower potential energy. 

Since the C20 fullerenes in the (020)2 dimer are, on the one hand, at small distances from 
each other and hence strongly interact with each other and, on the other hand, both preserve 
the overall spherical shape (see Fig. lb), the (C 2 o)2 dimer can loose its molecular structure 
both through the fusion of two C 2 o fullerenes to the C 40 cluster and through the decay of 
one of them [12]. A question arises of which of two factors, thermodynamic (the minimum of 
the energy of the final product after fusion) or kinetic (the thermally activated transition to 
the metastable configuration with a high energy upon decay), finally determines the stability 
of the (C 20 ) 2 dimer? To answer this question, we study the dynamics of the (C 2 o)2 dimer 
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numerically at high (sufficient to overcome the corresponding potential barriers U ~ 5 eV 
[12]) temperatures. 

To calculate the energies of arbitrary atomic configurations, we used the tight binding 
model [13], which is a reasonable compromise between ab initio calculations and too simplified 
empirical approaches with classical inter-atomic potentials. This model well describes both 
small clusters and macroscopic carbon structures. Previously, we used it to simulate the 
fullerenes C20, Ceo, and other carbon systems [10, 12, 14-18]. For the binding energy of an 
isolated C 20 fullerene, this method gives E b = 6.08 eV/atom [12] (in agreement with a value 
of E b = 6.36 eV/atom obtained in ab initio calculations [19). For the binding energy of the 
open-[2+2] isomer of the (C 2 o)2 dimer, see Fig. lb, this method gives a value of E b = 6.20 
eV/atom (which corresponds to the coagulation energy of two fullerenes AE = 2E pot [C 2 o\- 
Epot[(C2o)2] = 4.9 eV, to be compared with the value AE = 6.3 eV obtained by the density 
functional method [8]). 

To analyze the thermal (it can also be called kinetic) stability of the (C 2 o)2 dimer, we 
study its dynamics numerically for the initial temperature T ini = 2000 - 3500 K and fixed 
total energy E pot + Eki n = const, where E^ n is the kinetic energy in the center-of-mass syatem. 
In this formulation of the problem, the system temperature T is a measure of the energy of the 
relative motion of atoms (for more details, see [12, 17, 18]). The values of T were calculated 
by the formula (E kin ) = |/c B T(3A^ — 6), where kg is the Boltzmann constant, N = 40 is 
the number of atoms in the (C 2 o)2 dimer, and the angular brackets stand for averaging over 
10 3 — 10 4 steps of the molecular dynamics simulation (the time of one step is to = 2.72 • 10 16 
s). 

We observed two fundamentally different types of the T(t) dependence. They are shown 
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in Fig. 2. In the first case (Fig. 2a), the dimer temperature at t ps 4.4 ns drops stepwise 
by AT 200 K. The decrease in temperature points to a decrease in the kinetic energy 
and, hence, an increase in the potential energy (because E pot + E kin = const). Therefore, the 
(020)2 dimer passes to a less favorable configuration. Our analysis showed that this transition 
is caused by the decay of one of the C20 fullerenes, whereas the other fullerene preserved the 
cage shape, see Fig. 3a. The binding energy of the metastable state formed upon the decay 
is equal to (after the relaxation to the local minimum of the potential energy) E b = 6.14 
eV/atom, which is lower than the value of E b for the initial (C 2 o)2 dimer by AE b = 0.06 
eV/atom. This AE b value is consistent with the value AT ps 200 K, taking into account that, 
in thermal oscillations near the local minimum of the potential energy E™™ , the average 
kinetic energy of the system is approximately (up to anharmonic corrections) equal to the 
average increase in the potential energy AE pot = E pot — E™ t n , so that, since the total energy 
is constant, i.e., E pot + E kin = E™ n + AE pot + E kin ^ E™ n + 2E kin ^ E™ n + 3Nk B T, in 
the 1 — » 2 transition from the vicinity of one local minimum to the vicinity of the other local 
minimum, the relation E™ t n ' L + 3iVfc B Ti « E™ t n ' + SNk B T 2 is valid, from which it follows 
that AT = Ti - T 2 ^ AE b /3k B for £g n ' 2 > £g n>1 taking into account the definition of E b . 

On contrary, for the second case (Fig. 2b), the temperature at t ~ 4.4 ns increases stepwise 
by AT ps 200 K. This indicates an increase in the kinetic energy, i.e., a decrease in the 
potential energy. Therefore, the system passes to the energetically favorable configuration. 
However, since the initial open-[2+2] isomer has the minimum energy among all (C 2 o)2 dimers, 
this transition indicates that the C 2 o fullerenes loss their individuality and merge into the C 40 
cluster. We analized the atomic structure after the temperature jump and showed that this is 
indeed the case. The C 40 cluster formed upon fusion is shown in Fig. 3b. It can be considered 
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as one of nonclassical, strongly imperfect isomers of the C 40 fullerene: the C-C bonds in it 
form not only pentagons and hexagons, but also two 9-member rings. The binding energy of 
this cluster after relaxation, E b = 6.25 eV/atom, is higher than the value for the (020)2 dimer 
by AE b = 0.05 eV/atom, which is consistent with a value of AT pa 200 K upon fusion. 

We also observed the fusion of the C20 fullerenes into isomers with higher binding energies. 
Fig. 4 shows the T(t) dependence for the case, where an imperfect isomer of the C 40 fullerene 
with one heptagon and Ef, = 6.49 eV/atom is formed upon fusion at t pa 0.34 ns. At t pa 0.46 
ns, the annealing of the defect takes place and the C 40 fullerene with E b = 6.53 eV/atom is 
formed, on the surface of which the C-C bonds form 12 pentagons and 10 hexagons, see Fig. 
5. In this case, the temperature increases sequentially by AT pa 1100 and 150 K, which is 
consistent with changes AE b , see above. Transitions of the (C 2 o)2 dimer into other isomers 
of the C 40 fullerene were also observed. The general statistics is as follows (N k is the number 
of A;-gons with k > 6 in the C 40 isomer). Among 57 sets of the initial displacements and 
velocities of atoms that correspond to different initial temperatures T ini , the fusion into the 
C 40 cluster with A 9 = 2 occurs most frequently (16 times), see Fig. 3b. In addition, the 
clusters with N 10 = 2 and N 8 — 1 (six times); A 9 = N$ = N 7 = 1 (once); A 8 = 2 (twice); 
A 8 = A 7 = 1 (once); A 7 = 1 (once); and A 7 = 1 (once) are formed. In other 29 cases, the 
decay of one of the C 2 o fullerenes takes place (Fig. 3a). The further evolution of the system 
after it looses its molecular structure is as follows. The decay of one of the fullerenes in the 
(C 2 o)2 dimer is followed by the decay of the second fullerene, while fusion of two fullerenes into 
the C 40 cluster is followed by the decay of this cluster, so that the final product is always a 
quasi-two-dimensional or quasi-one-dimensional configuration as in the decay of the isolated 
C 2 o fullerene [12]. We never observed the return of the system to the initial open- [2+2] 
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configuration. The fission of the (C 2 o)2 dimer into two isolated C 2 o fullerenes has not been 
observed as well. 

There is correlation between the value of T ini and the mechanism of the loss of the stability 
of the (020)2 dimer. The loss of the molecular structure of the dimer occurs primarily due to 
fusion of two C20 fullerenes into the C40 cluster at T ini < 2500 K, due to the decay of one of 
the fullerenes at T ini > 3000 K, and due to both those processes at 2500 K < T ini < 3000 K, 
see Fig. 6. To determine the activation energies E a for the decay and fusion, we analyzed the 
temperature dependence of the lifetime r of the (C 2 o)2 dimer, using the Arrhenius formula 
with the finite-heat-bath correction [18, 20, 21] 

E a 



r -1 (T) = A-exp 



(1) 



k B T (1 - EJ2CT) 

Here, A is the frequency factor, C is the microcanonical specific heat, which we set to C — 
(37V - 6)k B , where N = 40 (as will be shown below, E a /2CT « 0.05 « 1 for T « 3000 
K). Figure 6 shows the calculated logarithm of r as a function of (T^w) -1 , where T* ni = 
Tiniil — E a /2CT ini ), see Eq. (1). It is seen that this dependence in a quite wide range 
T ini = 2000 — 3000 K is given, in the first approximation, by the straight line for both the 
decay and fusion processes; i.e., A(T) ~ const (since the spread of the points is sufficiently 
large, particularly for fusion processes, it is impossible to reliably determine the temperature 
dependence of A). Using the slopes of these straight lines and the points of their intersection 
with the ordinate axis, we determined the activation energies E a = 3.4 ± 0.2 eV and E a = 
2.7±0.3eV, as well as the frequency factors A = (1.8±0.3)-10 17 s" 1 and A = (5.6±1.1)-10 15 s" 1 
for the decay and coalescence processes, respectively. [Since E a enters both into the nominator 
of the exponential in Eq. (1) and into the renormalized initial temperature, it is determined 
self-consistently by successive iterations] . The values found for E a are one fourth to one third 



of the fragmentation activation energy of the C 60 fullerene (frequency factors are smaller by 
two to four orders of magnitude) [18]. 

In conclusion, we note that the value E a pa 3.4 eV for the decay of one C 20 fullerene in the 
(C 2 o)2 dimer is approximately one half of the value E a pa 7 eV for the decay of the isolated 
C20 fullerene [12]. This fact, along with the appearance of the second mechanism of the loss of 
the molecular structure (fusion of two C20 fullerenes into the C40 cluster) with a lower value of 
E a pa 2.7 eV, indicates that stability decreases when the (020)2 dimer is formed from two C20 
cages. In three-dimensional cluster structures based on C 2 o fullerenes, both the stabilization 
of the metastable state due to an increase in the number of nearest neighbors of each C20 
fullerene and the appearance of new channels of stability loss can be expected. We are going 
to study this problem in the next work. 
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Fig. 2. Time dependences of the temperature T of the (020)2 cluster dimer obtained from 
the molecular dynamics simulation by averaging the kinetic energy over each 10 4 steps. The 
total energy of the system is constant. The time of one step is equal to t = 2.72 ■ 1CT 16 s. 
The initial temperature is T ini = 2415 ± 10 K (a) and 2285 ± 10 K (b). 
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Fig. 3. (a) Atomic configuration after the decay of one of the C 20 fullerenes in the (C 2 o)2 
dimer, see Fig. 2a. The binding energy after relaxation is E b = 6.14 eV/atom. (b) The C 40 
cluster formed upon fusion of two C 2 o fullerenes in the (C 2 o)2 dimer, see Fig. 2b. The binding 
energy after relaxation is E b = 6.25 eV/atom. 
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Fig. 4. Same as in Fig. 2, for T ini = 2820 ± 10 K. 
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(a) 




Fig. 5. (a) Defect isomer of the C 40 fullerene that is formed upon fusion of two C 20 fullerenes 
in the (C 2 o)2 dimer, see Fig. 4. The binding energy after relaxation is E b = 6.49 eV/atom. 
(b) The C40 fullerene formed after annealing of the imperfect isomer. The binding energy 
after relaxation is E b = 6.53 eV/atom. The distant atoms are not shown for clarity. 
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Fig. 6. Logarithm of the lifetime r (in seconds) of the (020)2 cluster dimer for the processes 
of the decay of one of C 2 o fullerenes (circles) or fusion of two C 2 o fullerenes into the C 40 cluster 
(squares) versus the inverse temperature (Kelvins), taking the finite- heat-bath correcrtion into 
account, see text. The solid and dashed lines are the linear approximations for decay and 
fusion, respectively. 
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